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Abstract. We investigate an XY spin-glass model in which both spins and couplings 
evolve in time: the spins change rapidly according to Glauber- type rules, whereas 
the couplings evolve slowly with a dynamics involving spin correlations and Gaussian 
disorder. For large times the model can be solved using replica theory. In contrast to 
the XY-model with static disordered couplings, solving the present model requires two 
levels of replicas, one for the spins and one for the couplings. Relevant order parameters 
are defined and a phase diagram is obtained upon making the replica-symmetric 
Ansatz. The system exhibits two different spin-glass phases, with distinct de Almeida- 
Thouless lines, marking continuous replica-symmetry breaking: one describing freezing 
of the spins only, and one describing freezing of both spins and couplings. 

Recently various models with a coupled dynamics of fast Ising spins and slow 
couplings have been studied (see e.g. [|l], |5| and references therein). In addition to 
physical motivations, such as understanding the simultaneous learning and retrieval 
in recurrent neural networks or the influence of slow atomic diffusion processes in 
disordered magnetic systems, there is a more theoretical interest in such models in 
that they generate the replica formalism for a finite number of replicas n. Moreover, 
the replica number is found to have a physical meaning as the ratio of two temperatures 
(characterizing the stochasticity in the spin dynamics and the coupling dynamics, 
respectively). In this letter we extend the methods and results obtained for Ising spin 
models to a classical XY spin-glass with dynamic couplings, whose spin variables are 
physically more realistic than Ising ones. In addition, the XY model is closely related 
to neural network models of coupled oscillators, which provide a phenomenological 
description of neuronal firing synchronization in brain tissue. We solve our model upon 
making the replica-symmetric Ansatz, and calculate the de Almeida- Thouless (AT) lines 
(of which here there are two types), where continuous transitions occur to phases 
with broken replica symmetry. In doing so we also improve the calculations of M. As 
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in the Ising case we find two qualitatively different types of spin-glass phases. In one 
spin-glass phase the spins do freeze in random directions, but on the time-scales of the 
coupling dynamics these 'frozen directions' change. In the second spin-glass phase the 
spins as well as the couplings freeze, such that even on the large time-scales the 'frozen 
directions' of the spins remain stationary. 

We choose a system of iV classical two-component spin variables Si = (cos^, sin0j), 
% = 1...N, and symmetric exchange interactions (or couplings) J^, with a Glauber- 
type dynamics such that for stationary choices of the couplings the microscopic spin 
probability density would evolve towards a Boltzmann distribution, with the standard 
Hamiltonian H({Si}, {Jij}) = — J2k<e Jki.Sk- St and with inverse temperature (3 = T _1 . 
The couplings are taken to be of infinite range. They will now themselves be allowed 
to evolve in time in a stochastic manner, partially in response to the states of the spins 
and to externally imposed biases. However, we assume that the spin dynamics is very 
fast compared to that of the couplings, such that on the time-scales of the couplings, the 
spins are effectively in equilibrium (i.e. we take the adiabatic limit). For the dynamics of 
the couplings the following Langevin form is proposed (which is the natural adaptation 
to XY spins of the choices originally made in [[l], [3| for Ising spins): 

£ {St-S^+Ky _ . ggft) i<j = l N (1) 

The term (Si ■ Sj), representing local spin correlations associated with the coupling J^, 
is a thermodynamic average over the Boltzmann distribution of the spins, given the 
instantaneous couplings {Jke}- External biases = fiNBij serve to steer the weights 
to some preferred values (note: this notation follows that of Q). The are choosen 
to be quenched random variables, drawn independently from a Gaussian probability 
distribution with mean B /N and variance B/N. The decay term /xJy in (P is added 
to limit the magnitude of the couplings. Finally, the terms f]ij(t) represent Gaussian 
white noise contributions, of zero mean and covariance (Vij(t)Vki(t')) = 2T SikSji5(t — t'), 
with associated temperature T = f3~ l . Factors of iV are introduced in order to ensure 
non-trivial behaviour in the thermodynamic limit iV — > oo. 

The three independent global symmetries of our model, which can be expressed 
efficiently in terms of the Pauli spin matrices a x and a z , are the following: 

inversion of both spin axes : Si — » —Si for all i 

inversion of one spin axis : Si — > o~ z Si for all i (2) 

permutation of spin axes : Si — > cr x Si for all i . 

Upon using algebraic relations such as a x a z a x = —a z and o z a x a z = —a x we see that in 
the high T (ergodic) regime these three global symmetries generate the following local 
identities, respectively: 

(Si) = o, (Si ■ a x Sj) = 0, (S l -a z S J )=0. (3) 
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The equilibrium solution of the probability density associated with the stochastic 
equation (|l|) for the couplings follows from the fact that ([!]) is conservative, i.e. that it 
can be written as 

with the following effective Hamiltionian for the couplings: 

H({Jij}) = ~\ log MUij}) + l^N £ J 2 M BmJm • (5) 

P Z k<l k<£ 

In this expression Zp({Jij}) = T^,g .-, exp[pY^k<e JktSk • S(\ is the partition function 
of the XY spins with instantaneous couplings {«/#}. Thus the stationary probability 
density for the couplings is also of a Boltzmann form, with the Hamiltonian (|5|), and 
the thermodynamics of the slow system (the couplings) are generated by the partition 
function Zp = Jl\k<edJk£ exp[— f3H ({Jij})}, leading to (modulo irrelevant prefactors): 



Y[ dJu [^({J ii })] ?1 exp 

k<£ 



k<l 



(6) 



Finally, we define the disorder-averaged free energy per site / = —(fiN)" 1 (log Z^)b, in 
which (-)b is an average over the {B^}. In contrast to standard systems with frozen 
disorder, the replica number n is here given by the ratio n = (3//3, and can take any real 
non-negative value. The limit n — > corresponds to a situation in which the coupling 
dynamics is driven by the Gaussian white noise, rather than by spin correlations; in the 
limit n — > oo the influence of spin correlations dominates. 

We carry out the disorder average using the identity \ogZ^ = lim r ^ r _1 [^— 1], 



evaluating the latter by analytic continuation from integer r. Our system with partition 
function Z^ is thus replicated r times; we label each replica by a Roman index. Each of 
the r functions Zs, in turn, is given by (^), and involves Zp({Jij}) n which is replaced 
by the product of n further replicas, labeled by Greek indices. For non-integer n, 
again analytic continuation is made from integer n. Performing the disorder average 
in / results in an expression involving nr coupled replicas of the original system: 
{Si} — > {S® a }, with a = l...n and a = l...r. For iV — > oo this expression can 
be evaluated in the familiar fashion of replica mean-field theory [0, by saddle-point 
integration. This procedure induces the following order parameters: 
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The horizontal bar denotes thermal averaging over the coupling dynamics with fixed 
biases {-By}. Comparison with (|3|) shows that the order parameters m", and v^f 
measure the breaking of the global symmetries (|2|). For simplicity we choose Bq = 0. 
We make the usual assumption that, in the absence of global symmetry-breaking forces, 
phase transitions can lead to at most local violation of the identities Thus the 
latter will remain valid if averaged over all sites, at any temperature, which implies that 
m" = o and = v^ b 3 = 0. The spin-glass order parameters g^f , on the other hand, are 
not related to simple global symmetries, and serve to characterize the various phases. 

The final stage of the calculation is to make the replica symmetry (RS) Ansatz. Since 
observables with identical Roman indices refer to system copies with identical couplings, 
whereas observables with identical Roman indices and identical Greek indices refer to 
system copies with identical couplings and identical spins, in the present model the RS 
Ansatz takes the form g£f = 5 ab {6 a/3 + qi[l - 5 a p]} + q [l - 5 ab ] (note: S% ■ S% = 1). 
The remaining two order parameters are determined as the solutions of the following 
coupled saddle-point equations: 
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D . , Jdz P(z) [IojzZT- 1 h(z~) hjzxPE-^jBqo, 

q = dx P{x) \ j=L= } (7) 

JdzP(z)iI (z~)] n I (zxP~-iJlBq~) 



qi = J dx P{x) 



Jdz P(z) [Io(zZ)} [Ii(05)] 2 / O (^5-Vi^g« 



Jdz P(z) [I Q {z-Z)] n I Q {zx(3~-^J\Bq 



with the two short-hands J=1//j(3, E = (3J~(J+B)qi — \Bq G , with P(x)=xe~^ x 6[x\, 
and where the functions I n (x) are the Modified Bessel functions [Q. Their physical 
meaning is given by 

1 ^ , . 1 



i i 

It is clear that < qo < q\ < 1. 

We have studied the fixed-point equations (|7||D after having first eliminated the 
parameter redundancy by putting B = 1 and J = 3, which resulted in the phase 
diagram in the n-T plane as shown in figure |l|. In addition to a paramagnetic phase 
(P), where go = Qi = 0, one finds two distinct spin-glass phases: SGI, where qi > 
but q = 0, and SG2, where both q 1 > and q > 0. The SGI phase describes freezing 
of the spins on the fast time-scales only (where spin equilibration occurs); on the large 
time-scales, where coupling equilibration occurs, one finds that, due to the slow motion 
of the couplings, the frozen spin directions continually change. In the SG2 phase, on 
the other hand, both spins and couplings freeze, with the net result that even on the 
large time-scales the frozen spin directions are 'pinned'. The SG1-SG2 transition is 
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always second order. The transition SG1-P is second order for n < 2 (in which case its 
location is given by B + J = AT 2 ), but first order for n > 2. When n further increases 
to n > 3.5, the SGI phase disappears, and the system exhibits a first order transtion 
from P to SG2. 

Additional transitions occur, corresponding to a continuous breaking of replica 
symmetry. The stability of the RS solutions is, as always, expressed in terms of the 
eigenvalues of the matrix of second derivatives of quadratic fluctuations at the saddle- 
point ||. We have calculated all eigenvalues and their multiplicity following || []] (full 
details will be published elsewhere |§). We find two replicon eigenvalues: Xg, associated 
with the Greek replicas and Xr, associated with the Roman replicas. We have drawn the 
corresponding AT-lines, where Xg and Xr are zero, respectively, in the phase diagram 
as dashed-dotted lines. At low spin temperature T replica symmetry is broken with 
respect to the Greek replicas (Xg = 0), whereas at low coupling temperature T replica 
symmetry is broken with respect to the Roman replicas (Xr = 0). This is illustrated 
more clearly by drawing the phase diagram in the T-T plane, as in figure |[ This second 
figure also shows that there is no re-entrance from SGI to RSB, when T is varied for 
fixed T. 

Qualitatively the phase diagram of the present system is very similar to that of 
the Ising spin-glass with dynamic couplings as studied in ||, see figure |3], including 
the behaviour of both the Greek AT-line Xg = and the Roman AT-line Xr = 0. 
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n 

Figure 1. Phase diagram of the XY spin-glass with slow dynamic couplings, drawn 
in the n-T plane; for Bo = 0, B = 1 and J = 3. P: paramagnetic phase, q\ = qo = 0. 
SGI: first spin-glass phase, qi > and q — (freezing on spin time-scales only). SG2: 
second spin-glass phase, qi > and go > (freezing on all time-scales). AT lines: 
An = (Roman replicon), Xq = (Greek replicon). 
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Here our results differ from, and improve upon, those of || (which is why we present 
figure [|, rather than just refer to ||). The set of eigenvalues given in || turn out 
to satisfy only part of the relevant orthogonality conditions used in their calculation. 
The replica symmetric solution in SGI is always stable with respect to the Roman 
replicas. In fact, we can show analytically that the Roman AT-line coincides with the 
SG1-SG2 transition line. Our model, with dynamics on two different time-scales, is 
reminiscent of a simple XY model with one step replica-symmetry breaking (1RSB), 
and our eigenvalues formally resemble e.g. those describing the stability of the 1RSB 
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Figure 2. Phase diagram of the XY spin-glass with slow dynamic couplings, drawn 
in the T-T plane; for B = 0, B = 1 and J = 3. Further notation as in Figure |l|. 
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solution in the perceptron model of flOfl. 



In conclusion, we have solved a classical XY spin-glass model in which both the 
spins and their couplings evolve stochastically, according to coupled equations, but on 
different time-scales. The solution of our model in RS Ansatz is mathematically similar 
to that of the XY model with static couplings, but with one step RSB. Qualitatively, the 
phase diagram, which exhibits two different spin-glass types and both first and second 
order transitions, resembles closely that of the Ising spin-glass with dynamic couplings, 
provided appropriate adjustment of the calculation of the AT-line in || is carried out. 
Our calculation shows how the methods used for solving the Ising case can be easily 
adapted to more complicated spin types, and illustrates the robustness of the structure 
and peculiarities of phase diagrams describing the behaviour of large spin systems with 
dynamic couplings. 
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